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ABSTRACT 

We report on the results of a three-year program of coordinated X-ray and optical 
monitoring of the narrow-line Seyfert 1 galaxy NGC 4051. The rapid continuum 
variations observed in the X-ray spectra are not detected in the optical, although the 
time-averaged X-ray and optical continuum fluxes are well-correlated. Variations in 
the flux of the broad H/? line are found to lag behind the optical continuum variations 
by 6 days (with an uncertainty of 2-3 days), and combining this with the line width 
yields a virial mass estimate of ~ 1.1 x 10 6 Mq, at the very low end of the distribution 
of AGN masses measured by line reverberation. Strong variability of He n A4686 is also 
detected, and the response time measured is similar to that of H/3, but with a much 
larger uncertainty. The He n A4686 line is almost five times broader than H/3, and it is 
strongly blueward asymmetric, as are the high-ionization UV lines recorded in archival 
spectra of NGC 4051. The data are consistent with the Balmer lines arising in a low 
to moderate inclination disk-like configuration, and the high-ionization lines arising 
in an outflowing wind, of which we observe preferentially the near side. Previous 
observations of the narrow-line region morphology of this source suggest that the 
system is inclined by ~ 50°, and if this is applicable to the broad H/3-emitting region, 
a central mass of ~ 1.4 x 10 6 Mq can be inferred. During the third year of monitoring, 
both the X-ray continuum and the HeilA4686 line went into extremely low states, 
although the optical continuum and the H/3 broad line were both still present and 
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variable. We suggest that the inner part of the accretion disk may have gone into an 
advection-dominated state, yielding little radiation from the hotter inner disk. 

Subject headings: galaxies: active — galaxies: individual (NGC 4051) — galaxies: 
nuclei — galaxies: Seyfert X-rays: galaxies 

1. Introduction 

The broad components of permitted emission lines in the spectra of active galactic nuclei 
(AGNs) typically have velocity widths of a few to several thousands of kilometers per second. The 
defining characteristic of the subclass of AGNs known as narrow- line Seyfert 1 galaxies (NLSls) 
is that the broad components of their emission lines are much narrower (Vpwhm £ 2000 km s- 1 ) 
than is typical for broad-line objects (Osterbrock & Pogge 1985). NLSls are extreme AGNs in 
other respects as well — their UV-optical properties correlate well with the Boroson & Green 
(1992) primary spectral eigenvector identified in principal component analysis. In other words, 
NLS1 classification correlates well with strong optical Fen emission and weak [Om] AA4959, 5007 
emission (Boiler, Brandt & Fink 1996). While the possible importance of such extreme AGNs 
has been recognized for two decades (Davidson & Kinman 1978), interest in NLSls has increased 
recently as their unusual X-ray properties have come to light: they have unusually steep soft and 
hard X-ray spectra (Puchnarewicz et al. 1992; Boiler, Brandt, & Fink 1996; Brandt, Mathur, 
& Elvis 1997) and undergo rapid non- linear variability (Boiler et al. 1997). While rapid, large 
amplitude variability in the UV-optical has not been reported, existing data do not address well 
the relationship between the X-ray and long-wavelength variations. A good compilation of the 
observed properties of NLSls is given by Taniguchi, Murayama, & Nagao (1999). 

Possible explanations for the narrowness of the permitted lines include the following: 

1. NLSls have more distant-than-normal BLRs (Mason, Puchnarewicz, & Jones 1996; Wandel 
& Boiler 1998). The line widths are attributed to virial motion (M ~ V 2 r/G) in the vicinity 
of a supermassive black hole, but the orbital velocities are smaller than in more typical 
AGNs because of greater distance to the central source. 

2. NLSls are low-inclination (i.e., nearly face-on) systems (Osterbrock & Pogge 1985). In 
this model, the line widths are again due to orbital motion around the central black hole, 
and the bulk of the broad-line region (BLR) gas orbits in a common plane that is almost 
perpendicular to the line of sight, leading to relatively small Doppler widths. 

3. NLSls have relatively low black-hole masses, but high accretion rates. Again, the basic 
assumption is that the BLR motions are virial, but the central source has a lower mass. The 
luminosity can be kept relatively high by supposing that the accretion rate (relative to the 
Eddington rate) is correspondingly high in these sources. 
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The third of these explanations forms the current paradigm for NLSls, as it explains not only the 
narrow emission lines, but possibly also the steepness of the soft X-ray spectrum: the temperature 
structure of an optically thick, geometrically thin accretion disk is given by 

T(r) = 6.3 x 10 5 (M/M Edd ) 1/4 M~ 1/4 (r/iis)" 3 / 4 K, (1) 

where Mg is the black hole mass in units of 10 8 M Q and Rs is the Schwarzschild radius (Shakura 
& Sunyaev 1973). The temperature of the inner regions of the disk scales like M _1//4 , so the 
strength of the soft X-rays in NLSls might plausibly be ascribed to a low central mass and a high 
accretion-rate flow. 

One way to distinguish among various explanations of the NLS1 phenomenon is to measure 
the size of the BLR, which can be done via reverberation-mapping techniques (Blandford & 
McKee 1982; Peterson 1993; Netzer & Peterson 1997). It is usually assumed that the variability 
of an emission-line L(t) to a variable continuum C(t) can be linearized to the form 

SL(t) = f *(r) 5C(t-r) dr, (2) 



where VP(t) is the "transfer function," which depends on the geometry and reprocessing physics 
of the BLR. A representative time scale for response of a line can be found by a simple 
cross-correlation of the continuum and emission-line light curves. By convolving eq. (Q) with C(t), 
one obtains the cross-correlation function 

Fccf(t) = I ^(t')F A cf(t - t') dr', (3) 

where Facf( t ) is the continuum autocorrelation function (Penston 1991; Peterson 1993). The 
cross correlation lag r can be taken to be the light-travel time across the BLR, so the BLR size is 
given by r = ct. By combining this with the emission-line width Vf W hm> the mass of the central 
source can be inferred to be 

M = /F ™ hmCT , (4) 

where / is a factor of order unity that depends on the still unknown geometry and kinematics 
of the BLR. There is an implicit assumption that the gravitational force of the central object 
dominates the kinematics of the BLR; this is formally unproven, but at least in the case of the 
well-studied Seyfert 1 galaxy NGC 5548, the reverberation-mapping data are consistent with the 
required Vpwhm °c r~ 1 / 2 relationship (Peterson & Wandel 1999). Virial mass estimates based on 
reverberation-mapping data are now available for nearly 40 AGNs (Wandel, Peterson, & Malkan 
1999; Kaspi et al. 2000). 

Beginning in early 1996, we undertook a program of contemporaneous X-ray and optical 
spectroscopic monitoring of the galaxy NGC 4051, the only NLS1 galaxy in Seyfert's (1943) 
original list of high surface-brightness galaxies with strong emission lines. The X-ray variability 
characteristics of NGC 4051 are typical of the NLS1 class (Lawrence et al. 1987; McHardy et al. 
1995). 
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X-ray observations were made with the Rossi X-Ray Timing Explorer (RXTE), and optical 
spectra were obtained with ground-based telescopes, as described below. The purpose of this 
program has been twofold: 

1. To determine the nature of the relationship between the X-ray and UV-optical continuum 
variations. This is a particularly interesting question in the case of NGC 4051 since the 
X-ray flux dropped to an extremely low level in 1998 May (Uttley et al. 1999), towards the 
end of this campaign. 

2. To determine the BLR size and virial mass via reverberation techniques. 

In this contribution, we present the results of this program, and discuss their implications for the 
nature of the NLS1 phenomenon. 

2. Observations and Data Analysis 

2.1. X-Ray Observations 

The X-ray observations were made with the large area (0.7 m 2 ) Proportional Counter Array 
(PCA) on RXTE (Bradt, Rothschild & Swank 1993). The observations shown here, which are 
part of a continuing monitoring program, cover the period from 1996 March to 1998 June. The 
observations were scheduled initially to cover the largest range of variability time scales with the 
smallest number of observations in order to determine the X-ray power spectrum efficiently. The 
program therefore consisted of observations approximately every 7 to 10 days throughout the first 
year, followed by observations approximately every 2 weeks thereafter. In addition, during the 
first year, there were two more intensive monitoring periods: a two-week period of twice-daily 
observations and a four-week period of of daily observations. Each observation lasted for typically 
1 ksec. In 1996 December, there was also a period of 3 days during which NGC 4051 was observed 
for a total of 70 ksec. 

The PCA consists of 5 proportional counter units (PCUs) but typically only 3 PCUs (numbers 
0, 1 and 2) were operational and all count rates refer to the total counts from those 3 PCUs. 
Where other than 3 PCUs were in operation, the count rate has been normalized to 3 PCUs. 
The PCUs have 3 layers. Here we only include data from the upper layer as this layer provides 
the highest signal-to-noise ratio for photons in the energy range 2-20 keV where the flux from 
NGC 4051 is strongest. We used standard "good time interval" (GTI) criteria to select data with 
the lowest background. We reject data obtained when the Earth elevation angle was less than 
10°, when the pointing offset from the source was > 0.02°, or during passage through the South 
Atlantic Anomaly, or up to 5 minutes afterwards. The PCA is a non-imaging device with a field of 
view of FWHM ~ 1° and so the background which we subtract must be calculated from a model. 
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Here we use the FTOOL routine PCABACKEST V2.0c, with the new "L7" model, to calculate 
the background. 

The resultant 2-10 keV lightcurve is shown in the top panel of Fig. 1. Further details of the 
X-ray light curves and variability are given by Papadakis et al. (2000). For NGC 4051, 10 counts 
s _1 (2-10 keV), from 3 PCUs, corresponds to a flux of 4 x 10~ n ergs cm -2 s -1 . 



2.2. Optical Spectroscopy 

Optical spectroscopic observations were obtained between UT 1996 January 12 (Julian 
Date = JD2450095) and 1998 July 28 (JD2451022), covering three separate observing seasons. 
Observations were made with the Ohio State CCD Spectrograph on the 1.8-m Perkins Telescope 
of the Ohio State and Ohio Wesleyan Universities at the Lowell Observatory (data set A) and with 
the FAST spectrograph (Fabricant et al. 1998) on the 1.5-m Tillinghast Reflector of the Center 
for Astrophysics on Mt. Hopkins (data set B). A log of these observations is presented in Table 1. 
Column (1) gives the UT date of each observation and the Julian Date is given in column (2). 
The origin of the data (set A or B) appears in column (3). The projected size of the spectrograph 
entrance aperture was 5".0 x 7".5 (i.e., a slit width of 5".0 and a cross-dispersion extraction window 
of 7".5) for all set A spectra and 3".0 x 4".6 for all set B spectra. In each case, the slit was oriented 
in the east- west direction (i.e., the slit position angle was always 90°). The nominal spectral 
resolution was 9 A for all set A spectra and 5 A for all set B spectra. The wavelength coverage 
of each spectrum is given in column (4), and the file name is given in column (5). All of these 
spectra are publicly available on the International AGN Watch site on the World-Wide WebQ. 

The spectroscopic images were processed in standard fashion for CCD frames, including bias 
subtraction, dark-count correction when necessary, flat-field correction, wavelength calibration, 
and flux calibration based on standard-star observations. Since even under photometric conditions, 
AGN spectrophotometry is rarely more accurate than ~10%, the usual technique of flux calibration 
by comparison with standard stars is far too poor for AGN variability studies. We therefore base 
our flux calibration on a scale defined by the observed flux in the prominent narrow [O in] AA4959, 
5007 doublet. These lines originate in a low-density region that is more spatially extended than 
the BLR or the continuum source, which for practical purposes can be regarded as point sources. 
The larger light-travel time and long recombination time ensure that any narrow-line variations 
will occur only over much longer time scales than of interest in this experiment. We therefore 
assume that the [O in] lines are constant in flux, and use these to scale each spectrum. All spectra 
are scaled to a constant flux of F([0lll] A5007) = (3.91 ± 0.12) x 10~ 13 ergs s _1 cm~ 2 , which is 
based on the mean of ten spectra from data set A that were obtained under photometric observing 



lr The light curves and spectra are available at URL |http: / / www.astronomy.ohio-state.edu/~agnwatch / All publicly 
available International AGN Watch data can be accessed at this site, which also includes complete references to 
published AGN Watch papers. 
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conditions during the 1996 observing season. The [O in] A5007 flux measured from photometric 
spectra from subsequent years supports our assumption that this value can be assumed to be 
constant over the time scales of interest. All of the spectra are adjusted in flux to have this 
value of FQOm] A5007) by employing the spectral scaling software described by van Groningen 
& Wanders (1992). The process is as follows: spectra are adjusted in flux by a multiplicative 
constant that is determined by comparing each spectrum to a "reference spectrum" that has been 
formed by averaging all the highest-quality (i.e., typically signal-to- noise ratios S/N £ 30) spectra 
and scaling this mean spectrum to the adopted [O m] A5007 flux. All individual spectra are scaled 
relative to the reference spectrum in a least-squares fashion that minimizes the [O in] residuals 
in the difference spectrum produced by subtracting the reference spectrum from each individual 
spectrum. This program also corrects for small zero-point wavelength-calibration errors between 
the individual spectra, and takes resolution differences into account. 

At this point, measurements of each of the spectra are made. The continuum flux at ~ 5100 A 
(in the rest frame of NGC 4051, z = 0.002418, based on 21-cm emission [deVaucouleurs et al. 
1991]) is determined by averaging the flux in the 5090-5120 A bandpass (in the observed frame). 
The H/3 emission-line flux is measured by assuming a linear underlying continuum between 
~ 4772 A and ~ 5105 A, and integrating the flux above this continuum between 4820 A and 4910 A 
(all wavelengths in the observed frame). The long- wavelength cutoff of this integration band is 
chosen to avoid the Fen contamination underneath [Oiii] A4959. We also note that no attempt has 
been made to correct for contamination of the line measurement by the narrow-line component of 
H/3, which is of course expected to be constant. 

We have also measured the flux in the He n A4686 line. Only set B spectra are suitable for 
this measurement, since set A spectra do not extend shortward far enough to provide a suitable 
short-wavelength continuum point. The Hell flux was measured by adopting a linear underlying 
continuum between ~ 4447 A and ~ 4775 A, and integrating the flux above this continuum 
between 4613 A and 4772 A. 

We then compare the independent light curves from the two sets of data to identify small 
systematic flux differences between the sets, as we have done in many previous experiments (see 
Peterson et al. 1999 and references therein). We attribute these small relative flux offsets to 
aperture effects, although the procedure we use also corrects for other unidentified systematic 
differences between data sets. We define a point-source correction factor ip by the equation 

F(H/3) trU e = <^F(H/3) observed . (5) 

This factor accounts for the fact that different apertures result in different amounts of light loss for 
the point-spread function (which describes the surface-brightness distribution of both the broad 
lines and the AGN continuum source) and the partially extended narrow-line region (NLR). After 
correcting for aperture effects on the point-spread function to narrow-line ratio, another correction 
needs to be applied to adjust for the different amounts of starlight admitted by different apertures. 
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An extended source correction G is thus defined as 

observed 

G. (6) 

The value of G is essentially the nominal difference in the contaminating host-galaxy flux between 
the two spectrograph entrance apertures employed. 

This intercalibration procedure is accomplished by comparing pairs of nearly simultaneous 
observations from the two data sets to determine <p and G. In practice, the interval which we define 
as "nearly simultaneous" is two days or less, which means that in principle any real variability 
that occurs on time scales this short tends to be somewhat suppressed by the process that allows 
us to merge the two data sets. In this case, the adjustment has very little impact on the final 
results because nearly all of the pairs of data separated by two days or less are from data set B and 
thus have not been adjusted relative to one another. We find that the best-fit constants for set B 
relative to set A are tp = 0.982 ± 0.048 and G = (-1.243 ± 0.729) x 10" 15 ergs s" 1 cm" 2 A -1 . Both 
the H/3 and HeiiA4686 fluxes are adjusted as in eq. (|5|); the fact that the factor ip is so close to 
unity indicates that most of the NLR in this galaxy arises within a few arcseconds of the nucleus. 

The final continuum ^(5100 A) and H/3 emission-line fluxes are given in Table 2. Simultaneous 
(to within 0.1 day) measurements were averaged, weighted by the reciprocal of their variances. 

3. Analysis 

3.1. Continuum Variability 

The light curves listed in Table 2 are plotted in Fig. 1, along with contemporaneous 2-10 keV 
X-ray light curves. The data shown here span three observing seasons, beginning in 1996 January 
and ending in 1998 July. A summary of the general variability characteristics is given in Table 3. 
For the complete data base and for individual subsets of the data as given in column (1), columns 
(2) - (4) give respectively the number of individual observations and the average and median time 
intervals between them. The mean flux is given in column (5), and columns (6) and (7) give two 
widely used measures of the variability, F var , the root-mean-square (rms) fractional variability 
corrected for measurement error, as defined by Rodrfguez-Pascual et al. (1997), and i? m ax, the 
ratio of maximum to minimum flux, respectively. Both F vax and i? ma x parameters are affected by 
contamination of the measured quantities by constant-flux components; the optical continuum 
values are somewhat diluted by the constant contribution of the underlying host galaxy, and the 
emission-line values are affected by both narrow-line contributions and probably slowly varying 
Fen emission as well. In any case, these will have only a modest effect on F var and -R max and 
inspection of Table 3 shows clearly that the large-amplitude, rapid variations that characterize 
X-rays in NLSls are much less pronounced in the optical spectrum. 

While there is a clear lack of correlated short time-scale behavior of the X-ray and optical 
continua, the light curves in Fig. 1 suggest that a correlation on longer time scales is possible. To 
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test this quantitatively, we have suppressed the rapid variations by smoothing both the optical 
continuum and X-ray light curves with a rectangular function of width 30 days, as shown in Fig. 
3, similar to what was done by Maoz, Edelson, & Nandra (2000) in a comparison of X-ray and 
optical variability in the Seyfert 1 galaxy NGC 3516. Cross-correlation of the overlapping parts of 
these light curves using the methodology described in the next section yields a lag of the optical 
variations relative to the X-ray of r = 6l^ 2 days with a correlation coefficient r max = 0.74, i.e., 
the mean X-ray and optical fluxes are indeed correlated once the high-frequency variability is 
suppressed. The lag between variations in the two wavebands is highly uncertain, but consistent 
with zero or any small time lag expected in the continuum emitting region. 



3.2. Emission-Line Variability 

Comparison of the optical continuum and emission-line light curves shows that the variations 
in each are quite similar, indicating that the time delay between them is small. The time delay 
between continuum and emission-line variations can be quantified by cross-correlation of the light 
curves. During the first year, there is a period from JD2450183 to JD2450262 in which the light 
curves are well-sampled and the character of the variations permits an accurate cross-correlation 
measurement. In Table 3, we refer to this subset of data as "subset 1", and we plot this part of the 
light curves in an expanded form in Fig. 2. We cross-correlate the data shown here by using the 
interpolation cross-correlation function (ICCF) method of Gaskell & Sparke (1986) and Gaskell & 
Peterson (1987) and the discrete correlation function (DCF) method of Edelson & Krolik (1988), 
where in both cases, we employ the specific implementation described by White & Peterson (1994). 

The results of the cross-correlation analysis are summarized in Table 4 and the cross- 
correlation functions are shown in Fig. 4. For both H/3 and HeiiA4686, Table 4 gives the ICCF 
centroid T ccn t> and the location r pea k of the maximum value of the correlation coefficient r max . 
The centroid T ccn t is computed using all points near r pea k with values greater than 0.8r max - The 
uncertainties quoted for r cen t and r pea k are based on the model-independent Monte-Carlo method 
described by Peterson et al. (1998). 

By combining this lag with the Doppler width of the emission line, we can estimate a virial 
mass, as in eq. (4); for consistency with Wandel et al. (1999) and Kaspi et al. (2000), we use 
/ = 3/\/2 in eq. (4). Since the broad emission-line features are comprised of a number of different 
components (or contaminants), it is desirable to measure the Doppler width of only the variable 
part of the emission line. In order to isolate the variable part of the emission line and exclude 
constant components (such as contamination from the NLR), we measure the relevant line widths 
in the rms spectrum constructed from all the set B spectra in subset 1. The mean spectrum is 
constructed by averaging all N(= 18) set B spectra in subset 1, i.e. 

1 N 

* , W = JfY, F *W> (7) 
i=i 
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where F{(X) is the flux density (in ergs s 1 cm 2 A x ) at wavelength A in the ith spectrum. The 
rms spectrum is similarly constructed as 

1/2 

(8) 

The mean and rms spectra are shown as the top two panels in Fig. 5. The widths of the H/3 and 
Hen A4686 emission lines (full-width at half maximum, V FWH m) are given in Table 4, as Doppler 
widths in the rest frame of NGC 4051. A virial mass is then computed as described by Wandel, 
Peterson, & Malkan (1999). On the basis of the H/3 variations, a mass of l.ltoJ x 10 6 M is 
inferred; unfortunately however, the extremely large uncertainty in the He n A4686 lag renders the 
virial mass obtained from it not very enlightening, but it is consistent with the H/3 result. 

It is also important to keep in mind that because of the unknown geometry and kinematics 
of the BLR, the virial mass is reliable to only about an order of magnitude, i.e., the systematic 
uncertainties are much larger than the errors quoted here. However, an independent estimate 
of the inclination of the system has been made by Christopoulou et al. (1997) on the basis of 
the NLR morphology and kinematics. These authors model the NLR as an outflowing biconical 
region of inclination 50° and half-opening angle 23°. If the BLR and NLR axes are coaligned, then 
correcting the virial mass for inclination gives a central mass of 1.4tJ;g x 1O 6 M . 

There are a number of important features in the rms spectrum that deserve attention: first, 
the constant components, such as the [O in] AA4959, 5007 narrow lines that are so prominent in 
the mean spectrum, are absent, as expected, in the rms spectrum (except for weak residuals which 
reflect the accuracy to which accurate flux calibration can be achieved). Second, the weak broad 
wings of H/3 that can be seen in the mean spectrum are much weaker in the rms spectrum, i.e., the 
line core is more variable than the line wings. This could occur if the higher velocity material is 
much farther away from the ionizing source, or if some significant component of the high-velocity 
gas is optically thin (e.g., Shields, Ferland, & Peterson 1995); on physical grounds, we prefer the 
latter explanation. Third, the rms spectrum also shows that the optical Fe n emission varied little, 
if at all, during this period. The optical spectra of AGNs have broad blended Fen features that 
extend from just longward of H7 to just shortward of H/3, and over the range ~ 5100-5600 A. In 
the mean spectrum, these Fe 11 emission features are quite strong, which makes it hard to isolate 
the He 11 A4686 emission. The Fe 11 blends are, however, virtually absent in the rms spectrum. 
Fourth, He II A4686 is very prominent in the rms spectrum and is much broader than H/?. The 
rest-frame width of this line in the rms spectrum is ~ 5430 km s _1 , very typical of the line widths 
seen in normal Seyfert 1 galaxies. The centroid of the HeiiA4686 is strongly blueshifted relative to 
H/3, by about 1400 km s _1 . In order to demonstrate that these properties of the HeiiA4686 profile 
are real, we also computed the mean and rms spectra based on (a) set A spectra obtained at the 
same time and (b) set B spectra obtained during the 1997 monitoring season (Year 2). The rms 
spectra from these subsets, shown in the bottom two panels of Fig. 5, show the same He 11 A4686 
profile characteristics seen in the set B Subset 1 data shown in the second panel. 



a(A) = 



N — 1 



£ (F,(A) - F(A)r 
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4. Discussion 

4.1. The Continuum 

As described in the previous section, the X-ray and optical continuum variations are not 
closely coupled on short time scales. The X-ray continuum shows large scale variations on short 
time scales, as is typical of the NLS1 class. The rapid variations seen in the X-ray are not detected 
in the optical, as has previously been reported by Done et al. (1990) for this same source and by 
Young et al. (1999) for IRAS 13224—3809. However, if we average over the short time-scale flares, 
then the X-ray and optical continuum variations do seem to be coupled, though the time resolution 
of this experiment is insufficient to determine whether there is any lag between them at the level 
of days or less. The absence of strong coupling of the hard X-ray and optical continuum variations 
argues against reprocessing models in which hard X-rays are absorbed by a dense plasma (such as 
the accretion disk) and the energy is re-radiated at longer wavelengths (Guilbert & Rees 1988). 
The X-ray variations are more suggestive of localized flaring types of activity that may arise in a 
patchy corona above the accretion disk. 

It has already been pointed out based on these same RXTE data (Uttley et al. 1999) that 
the X-ray continuum of NGC 4051 virtually "turned off" in early 1998 (around JD 2450800; see 
Fig. 1). However, the optical spectroscopic data show that the optical continuum and emission 
lines (and therefore, by inference, the ionizing UV continuum) did not disappear at the same time. 
This seems in a sense rather suggestive of the kind of behavior that has been seen in Galactic 
black-hole systems such as GRO J 1655— 40 (Orosz et al. 1997). A possible interpretation of 
the behavior of NGC 4051 is that the inner X-ray producing (though not necessarily by purely 
thermal emission) part of the accretion disk in NGC 4051 has entered an advection-dominated 
accretion-flow (ADAF) state, in which radiation is emitted with very low efficiency (Narayan & 
Yi 1994; Narayan et al. 1998). This implies that there is a transition radius rt ra n inside of which 
the disk is an ADAF and outside of which it radiates efficiently, perhaps like a classical thin disk 
(Shakura &: Sunyaev 1973). The persistence of the optical continuum and emission lines suggests 
that this transition radius is somewhere between the regions that are most responsible for the soft 
X-rays and the H-ionizing continuum. We comment on this further in section 4.3 below. 

4.2. The Virial Mass and Implications for NLSls 

As noted above, reverberation-based size estimates for the broad emission lines and resulting 
virial mass estimates provide a potential means of distinguishing among the various NLS1 
models. In Fig. 6, we show the relationship between the BLR radius as measured from the 
H/? lag as a function of the optical continuum luminosity for all AGNs with Balmer line lags 
known to reasonable accuracy. All data used here are taken from the compilation of Kaspi et 
al. (2000), though their parameters for NGC 4051 are superceded by the values reported here. 
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This compilation contains six additional AGNs that could in principle be classified as NLSls as 
they meet the criterion 7 FWHM 2000 km s _1 . These sources, which we shall refer to below as 
"narrow-line objects", are the Seyfert galaxies Mrk 335 and Mrk 110 (from Wandel, Peterson, 
k Malkan 1999) and the QSOs PG 0026+120, PG 1211+143, PG 1351+640, and PG 1704+608 
(from Kaspi et al. 2000). 

The best-fit regression line (-Rblr °c ^o.62±o.02^ Dasec i on all objects from Table 5 of Kaspi et 
al. (2000) except NGC 4051, is shown as a dotted line. NGC 4051 lies approximately 2.8a above 
this regression line, although all the other narrow-line objects clearly fall in the locus defined by 
the AGNs with broader lines. It is difficult to argue that NGC 4051 is somehow different from the 
other AGNs, as there are several AGNs that have large displacements from the regression line. 
This is reflected in the large value of the reduced x 2 statistic, x 2 = 15.7, for this fit. 

In Fig. 7, we plot the mass-luminosity relationship for the AGNs from Kaspi et al., and 
we show (a) the best-fit regression line based on all objects except the seven narrow-line objects 
and (b) that based on the narrow-line objects alone. Formally, the slopes a (for M oc L a ) are 
significantly different, with a = 0.46 ± 0.06 for the narrow-line objects (including NGC 4051), and 
a = 0.27 ± 0.03 for the others. These two fits are separated by typically an order of magnitude 
in black-hole mass; the black holes in the narrow-line objects are about a factor of 10 lower than 
those of other AGNs of comparable luminosity. 

How well do these results allow us to distinguish among the various explanations for the NLS1 
phenomenon? We consider the possibilities: 

1. Do the BLRs of NLSls have anomalously large radii? The position of NGC 4051 in Fig. 6 
might suggest that this is possible, but the distribution of other narrow-line objects does 
not support this. Furthermore, as noted above, the scatter in the BLR-radius luminosity 
relationship is very large, and NGC 4051 is in a statistical sense not the largest outlier in 
this relationship (simply because other sources have smaller uncertainties in their measured 
lags). 

2. Are NLSls simply low-inclination systems? If the BLR is a flattened system, at low 
inclination (i.e., nearly face-on) the line widths will be decreased by a factor sini, but the 
measured emission-line lags will be relatively unaffected. On the other hand, assuming that 
the UV-optical continuum arises in an accretion disk at the same inclination, the apparent 
UV-optical luminosity is higher at lower inclination (e.g., see Fig. 32 of Netzer 1990). 
Thus, relative to similar sources at intermediate inclination, the masses of low-inclination 
sources will be underestimated, and their luminosities will be overestimated, displacing the 
narrow-line objects in Fig. 7 towards the lower right. This is generally consistent with the 
location of all seven of the narrow-line objects, including NGC 4051. The line transfer 
function for H/3 would provide a more definitive test of this hypothesis since it would allow 
determination of the inclination of the system. This would require more and better data 
than we have obtained in this experiment. 
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However, as noted earlier, Christopoulou et al. (1997) show that the NLR morphology and 
kinematics suggest a system that is inclined to the line of sight by ~ 50°. It seems reasonable 
to suppose that the NLR and BLR axes are approximately coaligned. If this is the case, 
then our virial estimate for the central mass is too small by a modest factor sin 50° ~ 0.77, 
and the corrected central mass is then 1.4+J;g x 10 6 M®. 

3. Are NLSls undermassive systems with relatively high accretion rates? Again, the distribution 
of the narrow-line objects, including NGC 4051, in Fig. 7 is consistent with this hypothesis. 
The narrow-line sources on this plot lie below the mass-luminosity relationship for other 
AGNs, at the lower end of the envelope around this relationship. 

In summary, the hypothesis that NLSls have unusually distant BLRs for their luminosity 
is probably not viable in general, although it could apply to the specific case of NGC 4051. At 
the present time, however, we cannot distinguish between the low-inclination and low-mass, high 
accretion-rate hypotheses on the basis of the reverberation results alone. The latter is favored 
on the basis X-ray considerations and the 50° inclination inferred from the NLR. Indeed, it is 
entirely possible that both effects (i.e., low inclination and low black-hole mass) contribute. An 
improvement in the optical spectroscopic monitoring data could allow determination of the H/3 
transfer function, which could allow discrimination between these competing models. 

4.3. Broad Hell Emission 

As noted in section 3.2, we detected very broad, blueshifted HeiiA4686 emission in the 
rms spectra. The blueshift of this component is suggestive of radial rather than virialized 
motion. There is no similar obvious component to the Balmer lines, but we should expect that 
similar blueshifted features might appear in the other higher-ionization lines in the UV. As no 
contemporaneous spectra were available, we retrieved the 31 International Ultraviolet Explorer 
SWP (Short- Wavelength Prime camera, wavelength range ~ 1150-2000 A) spectra from the Space 
Telescope Science Institute Multimission Data Archive. We formed an average of all these spectra, 
since the individual spectra were of rather low signal-to-noise ratio. In Fig. 8, we show the line 
profiles of Hen A4686 (the rms profile, as in Fig. 5) and those of Hen A1640 and Civ A1549 based 
on the mean IUE spectra. We note that each of these lines has a strong wing extending several 
thousand kilometers per second blueward of the line peak; indeed, the comparatively large widths 
and blueshifts of the UV lines in NLS1 spectra has been noted earlier by Rodriguez-Pascual, 
Mas-Hesse, & Santos-Lleo (1997). It is possible that this gas is related to the known warm 
absorber in NGC 4051 (McHardy et al. 1995). Interestingly, photoionization equilibrium modeling 
of the X-ray warm absorber data (Nicastro et al. 1999) suggests a distance from the source of 
approximately 5 light days, which is consistent with the reverberation result given in Table 4. 

The differences between the characteristics of the H/3 emission line on one hand and of the 
high ionization lines on the other suggests a two-component BLR, which has been proposed on 
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numerous occasions on other grounds (e.g., Collin-Souffrin et al. 1988; van Groningen 1987). In 
this particular case, an interpretation that is at least qualitatively consistent with all the data 
and relatively simple is that the Balmer lines arise primarily in material that is in a flattened 
disk-like configuration at a low to moderate inclination (to account for the narrow width of the H/3 
line), and the high-ionization lines arise in an outflowing wind, of which we see preferentially the 
component on the near side of the disk (to account for high velocity and blueward asymmetry). 
Such a model is illustrated schematically by Collin-Souffrin et al. (1988; their Fig. 1) and more 
recently by Dultzin-Hacyan, Taniguichi, & Uranga (1999; their Fig. 1). This geometry was also 
suggested by van Groningen (1987) to explain the line profiles and profile ratios of the Balmer 
lines in Mrk 335, another narrow-line object. 

As noted earlier, our virial mass estimate of M = 1.1 x 10 6 M Q might seriously underestimate 
the black hole mass if the inclination of the system is very low. Moreover, some low-inclination 
accretion-disk models predict relatively strong, variable EUV/soft X-ray fluxes (e.g., Netzer 1987; 
Madau 1988), consistent with observations of NLSls. We therefore cannot rule out the possibility 
that the NLS1 phenomenon is due principally to inclination effects. However, the strong rapid 
X-ray variability of NLSls seem to favor the low-mass, high accretion-rate explanation, as does 
the 50° inclination of the NLR, unless the BLR and NLR axes of symmetry are very different, 
which seems rather unlikely on physical grounds. 

The behavior of the He n A4686 line during Year 3 may provide additional information about 
the ionizing continuum at the time the X-rays went into an extremely faint state. In Fig. 9, we 
show the rms spectrum based on set B spectra obtained between JD2450810 and JD2451022. The 
H/3 line is strong and narrow, as it is in the other rms spectra shown in Fig. 5, which indicates that 
the continuum shortward of 912 A is still present and variable. However, the HeiiA4686 line is 
absent or very weak, indicating that the driving extreme ultraviolet (EUV) continuum shortward 
of the Hen edge at 228 A (54.4 eV) is not driving He n A4686 variations, either because the EUV 
flux is low or varying little. Evidence from earlier years (Uttley et al. 1999) shows that the 
EUV and X-ray fluxes in NGC 4051 are well correlated, which implies that the EUV continuum 
might also have been extremely weak in Year 3. Simultaneous BeppoSAX, RXTE, and EUVE 
observations obtained during the low state indicate that this is correct (Uttley et al. 2000). If the 
Hen emission has dramatically decreased during Year 3, we could infer that the transition radius 
between the inner ADAF region and the outer thin-disk structure must occur somewhere outside 
(or in the vicinity of) the region of the disk that contributes most of the flux at about 228 A, but 
inside the radius that contributes most of the flux at about 912 A. 

In this regard, whether or not there is residual He n A4686 emission in the low state during 
Year 3 becomes an interesting question. Note in particular that the measurements of the 
He II A4686 flux given in Table 2 and Figure 1 tell us little because these values include flux from 
broad-line Fe n blends and various narrow-line features. In order to address the question of how 
much He n A4686 persists in the low state, we have attempted a decomposition of the He n A4686 
feature. Our first step was to form average high-state (as in the top panel of Fig. 5) and low-state 
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spectra. We then attempted to remove the Fe n blends from the spectra by using the optical Fe II 
template spectrum of Boroson & Green (1982), convolved with a Gaussian to match the width of 
the features in the spectrum of NGC 4051. After subtraction of the flux-scaled and broadened 
Fe II template, we subtracted a power-law continuum from each spectrum. The He n A4686 region 
of the resulting high-state and low-state spectra are shown in the top panel of Fig. 10. We then 
subtracted the low-state spectrum from the high-state spectrum, forming the difference spectrum 
shown in the middle panel of Fig. 10. The flux in the Hen A4686 line in the difference spectrum is 
7.4 x 10 -15 ergs s _1 cm~ 2 . We then make the assumption that the Hen line profile in the difference 
spectrum can be used to model the He n contribution in the low-state spectrum. We then scale the 
He II profile in flux on a trial-and-error basis and subtract it from the low-state spectrum until the 
flux in the model profile causes the residual to have negative flux. Using this procedure, we find 
the largest possible value for the HeliA4686 line in the low-state is ~ 8.8 x 10~ 15 ergs s _1 cm~ 2 . 
These values mean that the He n A4686 flux decreased by at least 46% between the Year 1 high 
state and the Year 3 low state. Our conclusion is thus that while broad He n A4686 probably did 
not completely disappear during the low state in Year 3, the line flux did decrease dramatically 
and any Hen variations during this period were too small to detect. 

In any case, these results underscore the importance of multiwavelength observations of 
NLSls in very low states. Presumably comparison of the strength of the various UV emission 
lines in the low state relative to the values in the high state could lead to determination of r tra n- 
Unfortunately, no ultraviolet data are available during the present campaign. 

5. Summary 

On the basis of three years of combined X-ray and optical spectroscopic monitoring of the 
narrow-line Seyfert 1 galaxy NGC 4051, we reach the following conclusions: 

1. The rapid and strong X-ray variations that characterize narrow-line Seyfert 1 galaxies 
are detected in our X-ray observations of NGC 4051, but are not detected in the optical, 
consistent with previous findings. 

2. On time scales of many weeks and longer, there does appear to be a correlation between the 
X-ray and optical continuum fluxes. 

3. The variable part of the H/3 emission line has a Doppler width of ~ 1100 km s _1 , and a time 
delay relative to the continuum of about six days. Combining these quantities leads to a 
virial mass estimate of ~ 1.1 x 1O 6 M . If we assume that the inclination of the system is 
50°, as suggested by the NLR study of Christopoulous et al. (1997), then the mass of the 
central source is ~ 1.4 x 10 6 M Q . 

4. The He n A4686 emission line is strongly variable, although an accurate time delay cannot 
be measured. This line is about five times as broad as the H/3 line, and is strongly blueward 
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asymmetric, as are the UV high-ionization lines in this object. 

5. In the BLR radius-luminosity relationship, we find that narrow-line objects (those with 
V FWHM ~ 2000 km s _1 ) seem to fall on the same locus as AGNs with broad lines. 

6. In the virial mass-luminosity relationship, narrow-line objects populate the low-mass end of 
a rather broad envelope; they have virial masses typically an order of magnitude lower than 
other AGNs of similar luminosity. 

7. During the third year of this program, the hard X-ray flux decreased by approximately a 
factor of 10, and the broad-line Hen A4686 flux decrease by nearly a factor of two and did not 
show detectable variations during this low state. At the same time, the optical continuum 
and broad H/3 emission line were only slightly fainter than previously and continued to vary 
significantly. This suggests that the innermost part of the accretion disk went into an ADAF 
state, greatly reducing the production of high-energy continuum photons from the inner part 
of the accretion disk. 

A picture that is consistent with the emission-line characteristics is one in which the Balmer 
lines arise primarily in a disk-like configuration seen at low-to-moderate inclination and the 
high-ionization lines arise primarily in an outflowing wind. The high-ionization lines are blueward 
asymmetric because we see emission preferentially from the near side, with the far side at least 
partially obscured by the disk component (which might be an extension of the accretion disk). 
This geometry requires that the much of the high-ionization line flux arises in a region of scale 
comparable to the disk system that emits the low-ionization lines and that the disk system is at 
least partially opaque to the Hell line radiation. If NGC 4051 is typical of the NLS1 class, then 
it might be that NLSls are best described as low-mass, high accretion-rate systems, although 
the possible role of inclination cannot be discounted. Indeed, the full explanation of the NLS1 
phenomenon may involve both inclination and black-hole mass. 

For support of this work, we are grateful to the National Science Foundation (grant 
AST-9420080 to The Ohio State University). We thank the referee, M.-H. Ulrich, for suggestions 
that improved this paper. 
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Fig. 1. — Continuum and emission-line light curves for NGC 4051 between 1996 January and 1998 
June. The top panel shows the X-ray light curve for NGC 4051 in counts per second over the 
2-10 keV bandpass; one count per second corresponds to a flux of approximately 10~ 12 ergs cm~ 2 
s _1 over the 2-10 keV range. The second panel shows the optical continuum light curve in units 
of 10 -15 ergs s~ 1 cm -2 A -1 , as in Table 2. The third and fourth panels show, respectively, the H/3 
and Hen A4686 emission-line light curves, in units of 10~ 13 ergs s _1 cm -2 , also as in Table 2. 

Fig. 2. — The light curves in Fig. 1 are expanded to show the well-sampled region used in the 
optical continuum/ emission-line cross-correlation analysis ("subsetl") between JD 2450183 and 
JD 2450262. Units are as in Fig. 1. 

Fig. 3. — The mean X-ray (filled circles, flux scale on left axis) and optical continuum (open circles, 
flux scale on right axis) light curves from as in Fig. 1 smoothed with a rectangular function of width 
30 days. This shows when the rapid X-ray variability is suppressed by averaging, the X-ray and 
optical variations are correlated on time scales of several weeks and longer. 

Fig. 4. — Emission-line cross-correlation functions are shown, interpolation cross-correlation 
function as solid lines and discrete correlation function values as points with associated error 
bars. The top panel shows the 5100 A continuum-H/3 cross-correlation, along with the continuum 
autocorrelation function shown as a dotted line. The bottom panel shows the 5100 A continuum- 
He II A4686 cross-correlation. Peak and centroid values are given in Table 4. 

Fig. 5. — Mean (top panel) and root-mean-square (second panel) spectra of NGC 4051 based on 
individual spectra obtained during the well-sampled interval between JD 2450183 and JD 2450262 
("subset 1"), as in Fig. 2. The spectra are plotted in the observed frame, and in units of 10 -15 ergs 
s _1 cm _2 A _1 , as in Table 2. These are based only on the data from set B. The HeiiA4686 line 
is very broad (FWHM ~ 5430 km s _1 ) and blueshifted (by ~ 1400 km s _1 relative to the narrow 
emission lines); these characteristics are also seen in the rms spectrum based on contemporaneous 
spectra from set A (third panel) and on set B spectra obtained the next year (bottom panel). 

Fig. 6. — The relationship between the size of the broad Balmer-line emitting region and optical 
luminosity for AGNs. The filled circles are Seyfert galaxies from Wandel, Peterson, & Malkan 
(1999) and the open circles are QSOs from Kaspi et al. (2000). The large triangles those AGNs 
from the same sources whose broad lines have widths less than ~ 2000km s _1 , i.e., the Seyfert 
galaxies Mrk 335 and Mrk 110 and the QSOs PG 0026+120, PG 1211+143, PG 1351+640, and 
PG 1704+608. NGC 4051 is shown as a filled diamond. 

Fig. 7. — The relationship between the reverberation-based virial mass and and optical luminosity 
for AGNs. Date points are coded as in Fig. 6. 



-19- 



Fig. 8. — The top panel shows the profile of He n A4686 in the rms spectrum from Year 1 (1996), as 
in Fig. 5. The profiles of He n A1640 and C iv A1549 from mean IUE spectra are shown in the middle 
and bottom panels, respectively. Fluxes are in units of 10~ 15 ergs s _1 cm~ 2 A -1 , and in each panel 
the velocity scale is based on an assumed systemic redshift of z = 0.002418 (deVaucouleurs et al. 
1991). The blue wing appears more pronounced in the top panel because the rms spectrum enhances 
the variable part of the line. Also, the optical and UV observations are not contemporaneous. 



Fig. 9. — The rms spectrum based on the set B spectra obtained during Year 3 between JD2450810 
and JD2451022. Note that H/3 remains strongly variable, but in comparison to the rms spectra 
shown in Fig. 5, the HeiiA4686 line is much weaker, if present at all. 



Fig. 10. — The top panel shows the He n A4686 region of the mean high-state and low-state spectra, 
with Fe II emission and the underlying continuum subtracted. The middle panel shows the difference 
between the high-state and low-state spectra. The low-state spectrum is shown for reference (note 
in particular the strong He II A4686 narrow-line component in the low-state spectrum. The bottom 
panel shows our best estimate of the low-state spectrum with the largest possible contribution of 
broad Hen emission removed. Data are plotted in the rest frame of NGC 4051. The vertical axis 
is flux in units of 10 -15 ergs s _1 cm~ 2 A -1 . 
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Table 1. Log of Spectroscopic Observations 



UT 


Julian Date 




Range 


File 


Date 


(-2450000) 


Code 


(A) 


Name 


(1) 


(2) 


(3) 


(4) 


(5) 


1996 Jan 12 


95.0 


B 


3670- 


-7450 


n00095b 


1996 Jan 15 


98.0 


B 


3650- 


-7520 


n00098b 


1996 Jan 19 


102.0 


B 


3650- 


-7410 


n00102b 


1996 Jan 22 


105.1 


B 


3660- 


-7410 


n00105b 


1996 Jan 25 


108.0 


B 


3650- 


-7520 


n00108b 


1996 Jan 29 


112.0 


B 


3640- 


-7500 


n00112b 


1996 Feb 05 


118.8 


A 


4510- 


-5670 


n00118a 


1996 Feb 09 


122.8 


A 


4520- 


-5670 


n00122a 


1996 Feb 10 


124.0 


B 


3660- 


-7400 


n00123b 


1996 Feb 14 


127.8 


A 


4520- 


-5670 


n00127a 


1996 Feb 15 


129.0 


B 


3660- 


-7420 


n00129b 


1996 Feb 18 


132.0 


B 


3660- 


-7530 


n00132b 


1996 Feb 21 


134.9 


B 


3660- 


-7410 


n00134b 


1996 Feb 23 


136.8 


A 


4520- 


-5680 


n00136a 


1996 Mar 07 


149.8 


A 


4520- 


-5680 


n00149a 


1996 Mar 22 


164.8 


A 


4530- 


-5680 


n00164a 


1996 Mar 25 


167.7 


B 


3640- 


-7500 


n00167b 


1996 Mar 26 


168.8 


A 


4520- 


-5690 


n00168a 


1996 Mar 28 


170.8 


B 


3650- 


-7390 


n00170b 


1996 Apr 10 


183.6 


B 


3660- 


-7420 


n00183b 


1996 Apr 12 


185.6 


B 


3670- 


-7420 


n00185b 


1996 Apr 12 


185.7 


A 


4540- 


-5690 


n00185a 


1996 Apr 15 


188.6 


B 


3660- 


-7410 


n00188b 


1996 Apr 18 


191.6 


B 


3660- 


-7410 


n00191b 


1996 Apr 20 


193.6 


B 


3660- 


-7400 


n00193b 


1996 Apr 25 


198.7 


A 


4510- 


-5680 


n00198a 


1996 Apr 25 


198.8 


B 


3650- 


-7410 


n00198b 


1996 Apr 26 


199.6 


B 


3660- 


-7400 


n00199b 


1996 May 03 


206.7 


A 


4530- 


-5690 


n00206a 


1996 May 08 


211.6 


B 


3660- 


-7400 


n00211b 


1996 May 09 


212.8 


A 


4540- 


-5680 


n00212a 


1996 May 10 


213.6 


B 


3660- 


-7520 


n00213b 


1996 May 15 


218.6 


B 


3660- 


-7510 


n00218b 


1996 May 17 


220.8 


A 


4540- 


-5700 


n00220a 
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Table 1 — Continued 



UT 


Julian Date 




Range 


File 


Date 


(-2450000) 


Code 


(A) 


Name 


(!) 


(2) 


(3) 


(4) 


(5) 


IQQfi Mav 18 


221.7 


B 


3660- 


-7400 


n00221h 

11 VJ V J Z Z J_ YJ 


IQQfi Mav 22 


225.6 


B 


3670- 


-7430 


n00225h 

11 VJVJZZ o yj 


IQQfi Mav 24 


227 8 

ZZ I . O 


A 


4540- 


-5700 


n00227a 

11UUZZ 1 d 


IQQfi Mav 2fi 


22Q 7 


R 

1_J 


3670- 


-7430 


n0022Qh 

11UUZZ Zf YJ 


IQQfi Mav 30 


233 7 

ZOO. ( 


A 


4530- 


-5670 


n00233a 

11U \J ZjOOCXi 


IQQfi Tnn Dfi 

J-i/tyU O LI 11 UU 


240 6 


R 


3660- 


-7400 


n00240h 


IQQfi Tun f)7 

l^t^U *J LI 1 1 \J t 


241.7 


A 


4530- 


-5690 


n00241 a 


IQQfi Tnn 14 

-Li/i/U U LI 11 X 4: 


248 7 

Z'iO. 1 


A 


4540- 


-5690 


llVJUZ'iOCl 


IQQfi Tnn 1fi 


9^0 7 


R 

1 _) 


3600- 


-7540 




IQQfi Tnn 1Q 


9^3 fi 


R 

1 J 


3670- 


-7530 


nDD9 c i3L) 


IQQfi Tnn 1Q 


9^3 7 


A 


4520- 


-5680 


nDD9 c i3a 


IQQfi Tnn 93 


9^7 7 


R 

1 J 


3660- 


-7420 


nDD9 c i7L) 


IQQfi Tnn 95 

±Uc/\J O LI 11 ZjO 


95Q 6 

ZO .U 


R 


3670- 


-7550 


n0095Qh 

llUUiiUt/ YJ 


IQQfi Tnn 28 

lt7tyVJ O LI 1 1 ZjO 


2fi2 7 

ZL7Z . 1 


A 


4520- 


-5680 


n009fi9a 


IQQfi Tnl 17 


281 6 


B 


3670- 


-7520 


n00281b 


IQQfi Tnl 2D 

±C/C/\J O Lll ZU 


284 7 


R 


3670- 


-7530 


n00284h 


IQQfi Tnl 24 

A-tJiyXJ O Lll Z"-± 


288 6 


R 


3670- 


-7530 


n00288h 


1 QQfi Dpp 1 1 

lj jU 1 _y "L ±± 


42Q 


R 


3650- 


-7510 


n0042Qh 


1 QQfi Dor 1 7 


435 


R 


3660-7520 


n00435h 

llVJWiOO YJ 


1QQ7 Tan f)2 

U J 1 *J Clll UZ 


451 

'-to J- . VJ 


B 


3670- 


-7510 


n 0045 lb 


1QQ7 Tan DQ 

li/f 1 U Clll VJij 


458 


R 


3670- 


-7510 


n00458h 

IIUUIUO YJ 


1QQ7 Tan 1fi 

-Li7£/ ( U £111 ±U 


4fi5 


R 


4000- 


-7500 


n004fi5h 

llUU4:UtJ YJ 


1QQ7 Tan 31 

li/J 1 «J Clll O-L 


47Q Q 


A 


4720- 


-5990 


n0047Qa 


1QQ7 Tan 31 


480 


R 


4000- 


-7500 


n00480h 


1 nn7 T7, , i . no 
199/ teu 03 


480.1 


D 


4000- 


-7500 


n0048ob 


1997 Feb 06 


485.9 


B 


4000- 


-7530 


n00485b 


1997 Feb 09 


489.0 


B 


4000- 


-7500 


n00488b 


1997 Feb 14 


493.9 


A 


4380- 


-5550 


n00493a 


1997 Feb 14 


494.0 


B 


4000- 


-7520 


n00494b 


1997 Feb 27 


506.8 


A 


4360- 


-5500 


n00506a 


1997 Mar 02 


510.0 


B 


3660- 


-7520 


n00509b 


1997 Mar 12 


519.8 


B 


3660- 


-7500 


n00519b 


1997 Mar 13 


520.8 


A 


4290- 


-5820 


n00520a 


1997 Mar 14 


521.8 


B 


3670- 


-7310 


n00521b 
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Table 1 — Continued 



UT 


Julian Date 




Range 


File 


Date 


/ c\ A f r\r\r\r\\ 

(-2450000) 


Code 


(A) 


Name 


(!) 


( 2 ) 


( 3 ) 


(4) 


(5) 


1 997 Mar 2D 

U J 1 IV J. (XI _VJ 


527 8 


A 


4390- 


-5940 


n00527a 


1007 \ nr nq 

1 Jt/ I JJ1 VJ iJ 


547 6 


B 


3670- 


-7500 


n 00547b 


1 QQ7 Anr 19 


550 8 


R 

1 J 


3670- 


-7520 


n00550h 


1 QQ7 Anr H 


551 8 


R 


3680- 


-7540 


n00551h 


1QQ7 Anr 9Q 


567 6 


R 


3670- 


-7530 


n00567h 

11UL/(JU l YJ 


1 997 Mav D1 


569 6 

(JVJtJ . VJ 


B 


3670- 


-7540 


n 00569b 

11 VJ VJ O VJ t/ YJ 


iqq7 Mav D8 


576 6 

U 1 u.u 


B 


3650- 


-7500 


n00576h 

11VJVJO 1 VJ YJ 


1 QQ7 Mav 1 fl 


578 6 


R 

1 J 


3650- 


-7520 


n00578h 

11VJVJ(J 1 O YJ 


1 QQ7 Mav 1 9 


^80. 7 


R 


3660- 


-7520 


llUU(JOU YJ 


1 QQ7 Mav 1 4 


^89 fi, 


R 

1 J 


3650- 


-7530 


HUUtJOi/ YJ 


1 QQ7 Mav 9Q 


5Q7 7 


R 

1 J 


3660- 


-7520 


LlxJKJOcf i YJ 


1QQ7 Tnn D1 
±cfc/ 1 u liii \j ± 


600 6 


R 

1 J 


4000- 


-7530 


11UVJUVJU LJ 


1QQ7 Tnn 03 

±t?ty i O Llll VJ<J 


609 6 


R 

1 J 


3660- 


-7540 


n00609h 

11VJVJUVJZ YJ 


iqq7 T lln f)5 

1~'J*J I O LI 1 1 VJ(J 


604 6 

U V/ r . \) 


B 


3670- 


-7530 


n00604h 

11 VJ VJ VJ VJ t: U 


1QQ7 Tnn DQ 

±t?ty i U LI 11 Uf 


608 6 


B 

1_J 


3670- 


-7530 


n00608h 

11UVJUVJO YJ 


1 QQ7 Tnn 1 1 

Ii7t7 1 U Llll ±± 


61 6 


R 

1 J 


3670- 


-7540 


nDDfil Dh 

11UUU1UU 


1997 Tnn 98 

1 Jl/ l <J Llll 


627 6 

VJZj 1 . \) 


B 


3670- 


-7530 


n00697h 

11 VJ VJ VJ Zj 1 VJ 


1997 Tnl D1 

i~u*j i o lii y 1 1 


630 7 


B 


3670- 


-7520 


n00630h 

11VJVJVJ(JVJ YJ 


iqq7 t^i or 


635 6 

VJtJ<J . VJ 


B 


3670- 


-7530 


n00635h 

11VJVJVJ(J(J YJ 


iqq7 T^i 19 

1~U*J 1 O Lll 1a 


641 6 


B 


3670- 


-7530 


n00641h 

IIVJVJVJt; 1 YJ 


1997 Tnl 14 

1~U*J 1 O Lll Xt: 


643 6 


B 


3670- 


-7520 


n00643h 

IIVJVJVJ^O YJ 


iqq7 Nnv 99 

±i?tJ ( 1NLJV Z,_' 


775 


A 


4310- 


-5800 


n00775a 

11UU i i OcX> 


iqq7 Nov 94 


777 

f 1 f . L' 


R 

1 J 


3750- 


-7510 


n00777h 

11UVJ 1 1 I u 


iqq7 Nov 9q 


789 


R 


3750- 


-7500 


n00789h 


i nn7 t v,,, n/i 
iyy/ Dec 04 


7Q7 n 


T3 
D 


3750-7500 


nOO/8/b 


1997 Dec 27 


810.1 


B 


3670-7530 


n00810b 


1998 Jan 20 


834.0 


B 


3680- 


-7540 


n00834b 


1998 Jan 24 


838.1 


B 


3660- 


-7550 


n00838b 


1998 Jan 25 


839.0 


A 


4330- 


-5840 


n00839a 


1998 Jan 26 


840.1 


B 


3670- 


-7540 


n00840b 


1998 Jan 29 


843.0 


B 


3650- 


-7520 


n00843b 


1998 Feb 02 


847.0 


B 


3660- 


-7540 


n00847b 


1998 Feb 22 


867.0 


B 


3650- 


-7530 


n00867b 


1998 Feb 24 


869.0 


B 


3670- 


-7540 


n00869b 
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Table 1 — Continued 



T TT 

UT 


Julian Date 




Range 


tile 


TA j_ 

Date 


( o ,.\ v r\c\r\f\\ 

(—2450000) 


Code 


( A ) 


AT 


(i) 


(2) 


(3) 


( A \ 

(4) 


(5) 


1998 Feb 28 


873.0 


B 


3670-7530 


n00872b 


1998 Mar 03 


875.8 


A 


4300-5830 


n00876a 


1998 Mar 03 


876.0 


B 


3660-7530 


n00876b 


1998 Mar 04 


877.0 


B 


3660-7510 


n00877b 


1998 Mar 13 


885.8 


A 


4380-5910 


n00886a 

A L \J \J l_l KJ \J tAi 


1998 Mar 20 


892.8 


B 


3660-7510 


n00892b 


1998 Aor 03 


906.9 


B 


3630-7490 


n00906b 


1998 Aor 18 


921.9 


B 


3610-7470 


n00921b 


1998 Anr 23 


926.8 


A 


4320-5840 


n00927a 


1998 Anr 27 

1 Jl/U -ill J 1 Zj ( 


930.8 


B 


3620-7500 


n00930b 


1998 Mav 02 


935.9 


B 


3630-7490 


n00935b 


1998 Mav 03 

1 1/ / V.- 1 -.'Ill V V-/ ^ 


936.6 


B 


3620-7490 


n00936b 


1998 May 16 


949.6 


B 


3670-7510 


n00949b 


1998 May 28 


961.6 


B 


3660-7530 


n00961b 


1998 Jun 02 


966.7 


B 


3660-7540 


n00966b 


1998 Jun 16 


980.7 


A 


4310-5850 


n00981a 


1998 Jun 16 


980.7 


B 


3660-7540 


n00980b 


1998 Jun 19 


983.7 


B 


3650-7520 


n00983b 


1998 Jun 24 


988.6 


B 


3680-7530 


n00988b 


1998 Jun 27 


991.7 


B 


3660-7530 


n00991b 


1998 Jun 30 


994.6 


B 


3680-7520 


n00994b 


1998 Jul 15 


1009.6 


B 


3710-7510 


n01009b 


1998 Jul 18 


1012.6 


B 


3670-7530 


n01012b 


1998 Jul 25 


1019.6 


B 


3700-7500 


n01019b 


1998 Jul 28 


1022.6 


B 


3670-7530 


n01022b 



Codes for Data Origin 
A 1.8-m Perkins Telescope + Ohio State CCD Spectrograph 
B 1.5-m Tillinghast Reflector + FAST Spectrograph 



Table 2. Light Curves 



Julian Date F A (5100A) F(H/?) F(HeiiA4686) 

(-2450000) (lO^^ergss^^m" 2 !" 1 ) (10~ 13 ergs s" 1 cuT 2 ) 

(1) (2) (3) (4) 



95.0 


13.41 ±0.54 


4.17 ±0.17 


4.06 ± 0.45 


98.0 


13.88 ±0.56 


4.54 ±0.18 


4.22 ± 0.46 


102.0 


14.42 ± 0.58 


4.91 ±0.20 


4.53 ±0.50 


105.1 


13.46 ± 0.54 


4.54 ±0.18 




108.0 


13.75 ±0.55 


4.39 ±0.17 




112.0 


13.75 ±0.55 


4.42 ±0.18 


3.80 ± 0.42 


118.8 


14.01 ±0.28 


4.20 ±0.08 




122.8 


14.12 ±0.28 


4.67 ±0.09 




124.0 


13.74 ±0.55 


4.72 ±0.19 


3.81 ±0.42 


127.8 


13.51 ±0.27 


4.66 ± 0.09 




129.0 


14.37 ±0.57 


5.09 ±0.20 


3.97 ±0.44 


132.0 


14.91 ±0.60 


4.42 ±0.18 


3.88 ± 0.43 


134.9 


13.39 ±0.54 


4.63 ±0.19 


3.99 ± 0.44 


136.8 


14.13 ±0.28 


4.53 ±0.09 




149.8 


14.70 ± 0.29 


4.91 ±0.10 




164.8 


14.29 ±0.29 


4.71 ±0.09 




167.7 


14.66 ±0.59 


4.92 ±0.20 


3.86 ± 0.43 


168.8 


14.58 ±0.29 


4.75 ± 0.09 




170.8 


14.17 ±0.57 


5.06 ±0.20 


4.14 ±0.46 


183.6 


14.66 ±0.59 


4.84 ±0.19 


4.53 ± 0.50 


185.6 


14.51 ±0.58 


5.29 ±0.21 


3.86 ± 0.42 


185.7 


14.06 ± 0.28 


5.12 ±0.10 




188.6 


14.61 ±0.58 


5.17 ±0.21 


4.73 ± 0.52 


191.6 


14.63 ±0.58 


5.35 ±0.21 


4.78 ± 0.52 


193.6 


15.02 ±0.60 


5.75 ±0.23 




198.7 


14.43 ± 0.29 


5.63 ±0.11 




198.8 


13.54 ±0.54 


5.55 ±0.22 


4.47 ± 0.49 


199.6 


14.02 ±0.56 


5.54 ±0.22 


4.83 ±0.53 


206.7 


14.08 ± 0.28 


5.63 ±0.11 




211.6 


13.54 ±0.54 


5.48 ±0.22 


4.63 ±0.51 


212.8 


13.15 ±0.26 


5.36 ±0.11 




213.6 


13.16 ±0.53 


5.28 ±0.21 


4.73 ± 0.52 


218.6 


12.95 ±0.52 


5.09 ±0.20 


4.28 ± 0.47 


220.8 


12.27 ±0.25 


4.59 ± 0.09 





Table 2 — Continued 



Julian Date F A (5100A) F(H/3) F(HeiiA4686) 

(-2450000) (lO-^ergss-^m- 2 !- 1 ) (10~ 13 ergs s" 1 cm" 2 ) 

(1) (2) (3) (4) 



221.7 


13.34 ±0.53 


4.28 ±0.17 


3.08 ± 0.34 


225.6 


12.26 ±0.49 


4.37 ±0.17 


2.91 ±0.32 


227.8 


13.14 ±0.26 


4.14 ±0.08 




229.7 


13.91 ±0.56 


4.30 ±0.17 


3.77 ±0.41 


233.7 


12.84 ±0.26 


4.51 ±0.09 




240.6 


13.18 ±0.53 


4.88 ±0.19 


3.76 ± 0.41 


241.7 


13.29 ±0.27 


4.78 ±0.10 




248.7 


13.93 ±0.28 


5.64 ±0.11 




250.7 


12.99 ±0.52 


5.68 ±0.23 


4.59 ± 0.50 


253.6 


11.67 ±0.47 


5.22 ±0.21 


3.53 ±0.39 


253.7 


12.52 ±0.25 


4.94 ±0.10 




257.7 


12.09 ± 0.48 


4.47 ±0.18 


3.45 ± 0.38 


259.6 


11.92 ±0.48 


4.57 ±0.18 


3.35 ±0.37 


262.7 


12.41 ±0.25 


4.17 ±0.08 




281.6 


13.64 ±0.55 


5.25 ±0.21 


4.17 ±0.46 


284.7 


14.68 ± 0.59 


5.89 ± 0.24 




288.6 


13.25 ±0.53 


5.45 ±0.22 


4.74 ± 0.52 


429.0 


12.98 ±0.52 


5.84 ±0.23 


4.46 ± 0.49 


435.0 


13.53 ±0.54 


5.46 ± 0.22 


4.71 ±0.52 


451.0 


12.84 ±0.51 


5.51 ±0.22 


4.05 ± 0.45 


458.0 


13.25 ±0.53 


5.38 ±0.22 


4.34 ± 0.48 


465.0 


12.30 ±0.49 


4.74 ±0.19 


4.56 ±0.50 


479.9 


12.81 ± 0.26 


5.04 ±0.10 




480.0 


12.68 ±0.51 


5.10 ±0.20 


4.67 ±0.51 


483.1 


11.95 ±0.48 


4.69 ±0.19 


3.08 ± 0.34 


485.9 


12.41 ±0.50 


4.70 ±0.19 


3.22 ±0.35 


489.0 


14.03 ± 0.56 


4.92 ± 0.20 


3.99 ± 0.44 


493.9 


12.67 ±0.25 


4.99 ±0.10 




494.0 


13.53 ±0.54 


4.83 ±0.19 


3.30 ±0.36 


506.8 


11.46 ±0.23 


4.69 ± 0.09 




510.0 


12.30 ± 0.49 


4.37 ±0.17 


3.27 ±0.36 


519.8 


12.51 ±0.50 


4.74 ±0.19 


3.68 ± 0.41 


520.8 


12.50 ±0.25 


5.12 ±0.10 




521.8 


12.65 ±0.51 


4.68 ±0.19 


3.49 ± 0.38 
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Table 2 — Continued 



Julian Date F A (5100A) F(H/3) F(HeiiA4686) 

(-2450000) (lO^ergss^cm- 2 !- 1 ) (10~ 13 ergs s" 1 cm" 2 ) 

(1) (2) (3) (4) 



527.8 


10.92 ±0.22 


4.57 ±0.09 




547.6 


11.80 ±0.47 


4.52 ±0.18 


3.57 ±0.39 


550.8 


11.74 ±0.47 


4.47 ±0.18 


2.96 ±0.32 


551.8 


11.59 ±0.46 


4.42 ±0.18 


3.11 ±0.34 


567.6 


10.88 ± 0.44 


4.48 ±0.18 


3.76 ± 0.41 


569.6 


11.76 ±0.47 


4.31 ±0.17 


3.12 ±0.34 


576.6 


12.17 ±0.49 


4.12 ±0.17 


2.91 ±0.32 


578.6 


12.02 ± 0.48 


4.32 ±0.17 


3.24 ±0.36 


580.7 


11.97 ±0.48 


4.91 ±0.20 


3.28 ±0.36 


582.6 


12.10 ±0.48 


4.49 ±0.18 


3.35 ±0.37 


597.7 


12.57 ±0.50 


4.24 ±0.17 


2.92 ± 0.32 


600.6 


12.44 ±0.50 


4.15 ±0.17 


2.59 ±0.28 


602.6 


11.62 ±0.47 


4.61 ±0.18 


3.53 ±0.39 


604.6 


12.00 ±0.48 


4.50 ±0.18 


3.59 ± 0.40 


608.6 


11.96 ±0.48 


4.16 ±0.17 


3.39 ± 0.37 


610.6 


10.99 ± 0.44 


4.24 ±0.17 


2.80 ±0.31 


627.6 


11.26 ±0.45 


4.01 ±0.16 


2.92 ±0.32 


630.7 


11.14 ±0.45 


3.88 ±0.16 


2.49 ±0.27 


635.6 


11.87 ±0.47 


4.24 ±0.17 


2.37 ±0.26 


641.6 


12.33 ±0.49 


4.71 ±0.19 


3.20 ±0.35 


643.6 


10.75 ±0.43 


4.46 ±0.18 


2.83 ±0.31 


775.0 


13.29 ±0.27 


5.04 ±0.10 




777.0 


12.69 ±0.51 


5.18 ±0.21 


3.61 ±0.40 


782.0 


12.86 ±0.51 


5.41 ±0.22 


3.83 ± 0.42 


787.0 


13.49 ±0.54 


5.27 ±0.21 


4.20 ±0.46 


810.1 


14.06 ±0.56 


3.95 ±0.16 


2.97 ±0.33 


834.0 


12.00 ±0.48 


4.85 ±0.19 


2.87 ±0.31 


838.1 


11.44 ±0.46 


4.59 ±0.18 


2.99 ±0.33 


839.0 


11.89 ±0.24 


4.89 ±0.10 




840.1 


12.50 ±0.50 


4.91 ±0.20 


3.40 ±0.37 


843.0 


12.45 ±0.50 


4.87 ±0.19 


3.29 ±0.36 


847.0 


12.33 ±0.49 


4.82 ±0.19 


3.20 ±0.35 


867.0 


12.25 ±0.49 


4.20 ±0.17 


2.86 ±0.31 


869.0 


12.55 ±0.50 


3.85 ±0.15 


2.53 ±0.28 
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Table 2 — Continued 



Julian Date 


Fa(5100 A) 


F(W) 


F(He II A4686) 


(-2450000) 


(10 10 ergs s 1 cm Z A 1 ) 


(10 13 ergs s 1 cm 2 ) 


(i) 


(2) 


(3) 


(4) 


873.0 


12 15 ± 49 


4.09 ±0.16 


3.01 ±0.33 


875.8 


12 29 ± 25 


4.49 ± 0.09 




876.0 


12 58 ± 50 


4.31 ±0.17 


3.09 ±0.34 


877.0 


19 44 ± 50 


4.37 ±0.17 


3.09 ±0.34 


885 8 

COO . o 


11 54 + n 23 


4.18 ±0.08 




8Q2 8 


i o 1 4 + n 4Q 


4.54 ±0.18 


3.05 ± 0.34 


906.9 


11 44 + 46 


4.53 ±0.18 


2.60 ±0.28 


921.9 


11 61 ± 46 


4.53 ± 0.18 


2.93 ±0.32 


926.8 


11 13 ± 22 


4.43 ±0.09 




930.8 


11.43 ± 0.46 


4.01 ±0.16 


2.74 ± 0.30 


935.9 


11.92 ±0.48 


4.51 ±0.18 




936.6 


11.23 ±0.45 


4.39 ±0.18 


2.65 ±0.29 


949.6 


12.03 ± 0.48 


5.04 ±0.20 


3.41 ± 0.38 


961.6 


11.32 ±0.45 


4.90 ±0.20 


3.09 ±0.34 


966.7 


11.60 ±0.46 


4.39 ±0.18 


3.00 ± 0.33 


980.7 


11.72 ±0.21 


4.73 ±0.08 


3.13 ±0.34 


983.7 


11.41 ±0.46 


4.46 ±0.18 


2.60 ±0.29 


988.6 


11.56 ±0.46 


4.30 ±0.17 


2.78 ±0.31 


991.7 


11.35 ±0.45 


4.15 ±0.17 


3.08 ± 0.34 


994.6 


11.37 ±0.46 


4.11 ±0.16 


2.56 ±0.28 


1009.6 


10.98 ±0.44 


3.76 ±0.15 


2.50 ±0.28 


1012.6 


12.18 ±0.49 


3.70 ±0.15 


2.98 ± 0.33 


1019.6 


12.75 ±0.51 


4.20 ±0.17 


3.33 ±0.37 


1022.6 


12.09 ± 0.48 


4.57 ±0.18 


3.02 ± 0.33 
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Table 3. Variability Parameters 



Subset 


No. Epochs 


At ave (days) 


Aimed (days) 


Mean flux a 


± var 


^max 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


X-Ray (2-10 kev): 












All 


140 


7.5 


7.3 


6.20 ±3.85 


0.618 


42.85 ± 16.36 


Year 1 


33 


2.8 


0.7 


7.87 ±3.25 


0.411 


5.766 ± 0.459 


Year 2 


20 


11.5 


13.0 


6.42 ±4.95 


0.770 


17.68 ± 2.03 


Year 3 


19 


14.1 


14.1 


1.97 ± 1.81 


0.906 


14.79 ± 5.65 


C!nntinnnm 


(5100 AV 












All 


126 


7.4 


3.1 


12.74 ± 1.08 


0.077 


1.397 ±0.079 


Year 1 


01 


3.9 


3.0 


13.66 ± 0.82 


0.050 


1.287 ±0.073 


Year 2 


38 


5.8 


3.2 


12.16 ±0.76 


0.050 


1.305 ± 0.074 


Year 3 


37 


6.9 


4.5 


12.06 ± 0.69 


0.043 


1.281 ±0.072 


Subset 1 


29 


2.8 


2.2 


13.38 ±0.92 


0.059 


1.287 ±0.073 


H/J: 














All 


126 


7.4 


3.1 


4.77 ± 0.48 


0.095 


1.592 ±0.092 


Year 1 


51 


3.9 


3.0 


4.91 ± 0.48 


0.092 


1.423 ±0.064 


Year 2 


38 


5.8 


3.2 


4.65 ± 0.44 


0.086 


1.505 ± 0.086 


Year 3 


37 


6.9 


4.5 


4.50 ± 0.42 


0.086 


1.462 ± 0.084 


Subset 1 


29 


2.8 


2.2 


5.02 ±0.51 


0.096 


1.389 ±0.062 


Hen A4686 














All 


93 


10.1 


4.5 


3.50 ±0.67 


0.154 


2.038 ±0.316 


Year 1 


29 


6.9 


4.0 


4.08 ±0.52 


0.064 


1.660 ±0.258 


Year 2 


33 


6.7 


3.6 


3.42 ± 0.62 


0.143 


1.987 ±0.309 


Year 3 


31 


8.2 


5.0 


3.04 ± 0.38 


0.057 


1.680 ± 0.263 


Subset 1 


17 


4.8 


3.5 


4.08 ± 0.65 


0.114 


1.660 ± 0.258 



a X-ray fluxes in counts per second in 2-10 keV band. Other fluxes in units as in Table 2. 
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Table 4. Cross-Correlation Results and Virial Mass Estimates 



Parameter 


F(W) 


F(Hell A4686) 


(1) 


(2) 


(3) 



Tcent (days) 


5.92±f;i| 


4 - 4y -5.60 


rpeak (days) 


4-6±f;i 


^ q+8.9 


'"max 


0.800 


0.767 


Vfwhm (kms" 1 ) 


1110 ± 190 


5430 ± 510 


M vir (10 6 M ) 


i 1+0.8 
i - i -0.5 


19.41^ 



I 

, ■ 

• -A 

' 'i.i 
I i • 1 1 1 ' ■ ■ i , ' • 




i i i i i i i i i i i i i i i i i i 

200 400 600 800 1000 
Julian Date - 2450000 



■ i 




20 
10 


16 
14 
12 
10 





• 
• 

•• 

• • * •%• V . 

••I * 
% f • 


• • 

• — 


1 1 


1 1 1 1 


1 


t 






- i,h i 

1 1 


1 1 1 1 


1 



200 250 
Julian Date - 2450000 




200 400 600 800 



Julian Date - 2450000 




_l I I I I I I I I I I 

-20 20 

Lag (days) 




I I I I I I I I I I I I I I H 

4600 4800 5000 5200 

o 

Wavelength (A) 



1000 F — | I I I I 1 1 1 1 1 I I I I 1 1 1 1 1 I I I I 1 1 1 1 1 I I I I 1 1 1 1 II I I I I btj 



100 = 



"O 



CD 



10 = 



DC 
_l 
CD 



DC 



1 = 





tt m\ 




•"0 



A ° ~ 







r ^ 




4 



f o 



MM 



I I I I 



10 41 10 42 10 43 10 44 10 45 10 46 

-1 



XL (5100 A) (ergs s"') 




XL X (5100 A) (ergs s" 1 ) 



.2 




-8000 -6000 -4000 -2000 2000 4000 

Velocity (km/sec) 




4600 4700 4800 

o 

Wavelength (A) 



